The Isotropic Radio Background and Annihilating Dark Matter 



(N 

o 

(N 



in 



O 
U 

6 

u 



> 

in 

en 

en 

o 

(N 

> 

X 



Dan Hooper^'^, Alexander V. Belikov^, Tesla E. Jeltema**, Tim Linden"', Stefano Profumo^, and Tracy R. Slatyer^ 

^ Center for Particle Astrophysics, Fermi National Accelerator Laboratory, Batavia, IL 60510, USA 

^Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA 

^Institut d'Astrophysique de Paris, UMR 7095, CNRS, 

UPMC Univ. Paris 06, 98 bis boulevard Arago, 75014 Pans, France 

"^Department of Physics and Santa Cruz Institute for Particle Physics, 

University of California, 1156 High Street, Santa Cruz, CA 95064, USA and 

'^School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, USA 

(Dated: February 27, 2013) 

Observations by ARCADE-2 and other telescopes sensitive to low frequency radiation have re- 
vealed the presence of an isotropic radio background with a hard spectral index. The intensity 
of this observed background is found to exceed the flux predicted from astrophysical sources by a 
factor of approximately 5-6. In this article, we consider the possibility that annihilating dark mat- 
ter particles provide the primary contribution to the observed isotropic radio background through 
the emission of synchrotron radiation from electron and positron annihilation products. For rea- 
sonable estimates of the magnetic fields present in clusters and galaxies, we find that dark matter 
could potentially account for the observed radio excess, but only if it annihilates mostly to electrons 
and/or muons, and only if it possesses a mass in the range of approximately 5-50 GeV. For such 
models, the annihilation cross section required to normalize the synchrotron signal to the observed 
excess is av ~ (0.4 — 30) x 10~^® cm^/s, similar to the value predicted for a simple thermal relic 



{av 



3 X 10" 



cm /s). We find that in any scenario in which dark matter annihilations are 



responsible for the observed excess radio emission, a significant fraction of the isotropic gamma ray 
background observed by Fermi must result from dark matter as well. 

PACS numbers: 95.85.Bh, 95.85.Fm, 95.35. -l-d; FERMILAB-PUB-12-072-A 



I. INTRODUCTION 

Over the past several years, evidence indicating the 
existence of a significant isotropic radio background has 
been uncovered by a number of instruments. In 2009, the 
ARCADE 2 (Absolute Radiometer for Cosmology, As- 
trophysics and Diffuse Emission) collaboration reported 
measurements of the absolute sky temperature at a num- 
ber of frequencies between 3 and 90 GHz [1, . While these 
measurements are dominated by the cosmic microwave 
background (CMB) at frequencies above several GHz, 
they reveal the presence of significant excess power at 
the lowest measured frequencies [3j. This conclusion is 
further supported by a number of observations at lower 
frequencies, as reported by Roger et al. (0.022 GHz) 3J, 
Maeda et al. (0.045 GHz) [4], Haslam et al. (0.408 
GHz) [g, and Reich and Reich (1.42 GHz) [B]. The emis- 
sion observed by each of these groups is in significant 
excess of what can be attributed to Galactic emission, 
or to unresolved members of known extragalactic radio 
source populations [U [5] . 

In the left frame of Fig. [l] we plot the measured spec- 
trum of the isotropic radio background (after removing 
galactic foregrounds and resolved radio sources) [5]. At 
frequencies below a few GHz, these measurements sig- 
nificantly exceed the uniform temperature of the CMB 
(which is shown as a horizontal dotted line). In the 
center frame, this spectrum is shown after subtracting 
the contribution of the CMB, revealing a power-law like 
spectrum with an index of approximately —2.6 (—1.6 
in dNj/dE.y units) extending over at least three orders 



of magnitude in frequency, between ~20 MHz and ^10 
GHz [T . And while a fraction of this emission is expected 
to originate from faint and currently unresolved radio 
sources, estimates of this contribution based on deep sur- 
veys from the Very Large Array (VLA), among other 
observations, find that such sources should account for 
only approximately 20% of the observed isotropic radio 
emission. In the right frame of Fig. [l] we show the spec- 
trum after subtracting the estimated contribution from 
unresolved radio sources, as described in Ref. [5]. 

The origin of the excess isotropic radio emission is cur- 
rently unknown. Various systematic effects to explain 
this excess have been ruled out [1], and galactic origins 
(such as free-free or galactic synchrotron emission) are 
difficult to accommodate [H [71 |H] . Standard astrophysi- 
cal sources, such as radio quiet quasars, radio supernovae, 
and diffuse emission from clusters, have been considered 
but appear to be unable to account for the excess emis- 
sion [Hiin]- 

Any astrophysical objects capable of producing the ob- 
served isotropic excess without exceeding the observed 
number of resolved radio sources must be both highly nu- 
merous and very faint [SJ [TU] . The best known example 
of such a source class are star-forming galaxies. How- 
ever, in order for star-forming galaxies to account for the 
entire isotropic radio flux without exceeding the amount 
of emission observed at infrared (IR) wavelengths, the 
ratio of radio-to-IR emission must be increased by ap- 
proximately a factor of 5 at high redshifts over what 
is observed from local sources. Such strong evolution 
seems in conflict with current measurements [111 . Fur- 
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FIG. 1: Measurements of the isotropic radio background, as reported by Roger et al. (0.022 GHz) [3j, Maeda et al. (0.045 
GHz) iJ, Haslam et al. (0.408 GHz) 5 , Reich and Reich (1.42 GHz) p;, and by the ARCADE 2 coilaboration (3.20 to 90 
GHz) [T] (see also, Ref. 12]). In the left frame, the total observed background is shown, including the contribution of the CMB 
(shown as a dotted line). In the middle frame, the CMB has been subtracted, revealing a power-law spectrum of isotropic 
emission. In the right frame, an estimate of the contribution from unresolved radio sources has been subtracted as well. See 
text for more details. 



theriTLore, measurements of the isotropic gamma ray flux 
also strongly constrain such a scenario 



An alternative possibility for the origin of the excess 
isotropic radio emission is annihilating dark matter [131 
114] . Electrons and positrons produced in dark matter an- 
nihilations lose energy through both synchrotron and in- 
verse Compton processes, the former of which contributes 
to the radio background. Radio synchrotron emission has 
long been considered as a probe of dark matter annihi- 
lations in a number of contexts: from the Inner Milky 
Way [16 , the Galactic Center [T7], the Milky Way's non- 
thermal radio filaments [TS], galaxy clusters [12], and 
nearby dwarf galaxies and other substructures [SD] . Dark 
matter halos are predicted to be very numerous and faint 
radio sources - every dark matter halo large enough to 
contain a significant magnetic field constitutes such a 
source - making them ideal candidates to generate the 
observed radio background [l^. And while individual 
extragalactic halos are not expected to be bright enough 
to be resolved as radio sources, the sum of all halos 
could collectively constitute a bright source of isotropic 
radio emission. Furthermore, unlike most astrophysical 
sources of radio emission, some models of dark matter 
annihilation produce comparatively little emission at IR 
or gamma-ray wavelengths, and thus are relatively un- 
constrained by such observations (especially in the case 
of annihilations to e^e^ or fi^fi^). 

In this paper, we calculate the spectrum of radio emis- 
sion from dark matter annihilations taking place within 
the halos of galaxies and galaxy clusters, and determine 
in what cases such annihilations could potentially ac- 
count for the observed isotropic radio excess. We find 
that among a relatively narrow range of dark matter 
models, it may be possible to explain the observed emis- 
sion. In particular, to produce the observed radio spec- 
trum while also evading constraints from gamma-ray ob- 
servations (of both the isotropic gamma ray background 
and local dwarf galaxies) , we find that we must consider 



models in which the dark matter particles are somewhat 
light (rriDM ^4 — 50 GeV) and annihilate significantly 
to either e^e~ or ^^ n~ final states. In such models, the 
overall normalization of the radio excess requires dark 
matter annihilation cross sections which are typically 
similar to the value predicted for a simple thermal relic 
[av ^ 3 X 10~^^ cm^/s). In all of the dark matter models 
we found to be capable of producing the observed radio 
excess, dark matter annihilation products are also pre- 
dicted to make up a significant fraction of the isotropic 
gamma ray background. 



II. THE CONTRIBUTION FROM DARK 

MATTER ANNIHILATIONS TO THE 

EXTRAGALACTIC RADIO BACKGROUND 

In most models of weakly interacting massive par- 
ticles (WIMPs), dark matter annihilations deposit en- 
ergy into a combination of Standard Model particles, in- 
cluding gamma rays, neutrinos, protons, electrons, and 
their antiparticles. Energetic electrons and positrons 
produced through this process subsequently lose energy 
through a combination of synchrotron and inverse Comp- 
ton scattering. While inverse Compton scattering leads 
to the production of gamma rays and X-rays [21] [22], 
synchrotron processes with micro-Gauss-scale magnetic 
fields yield photons at radio and microwave wavelengths. 

The dark matter annihilation rate per volume, at red- 
shift z, is given by: 



i^W-/^(M,.)(l + .)3dM^ 



DM 



(1) 

where dn/dM is the differential comoving number den- 
sity of dark matter halos of mass M, av is the dark mat- 
ter annihilation cross section (multiplied by the relative 
velocity of the two WIMPs), mDM is the mass of the dark 



matter particle, and p(x, M, z) is the density of dark mat- 
ter at a location x within a halo of mass M. 

We calculate the number density of dark matter ha- 
los of a given mass as a function of redshift, j|j(Af, z) 
as described in Ref. [12], using standard values for cos- 
mological parameters. To describe the distribution of 
dark matter within a given halo, we adopt a generalized 
Navarro-Frenk- White (NFW) form [53] (7 = 1 being the 
canonical NFW value): 
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where Rg is the scale radius of the halo, which is related 
to the virial radius in terms of the halo concentration, 
c = i?vir/^s- For default parameters, we adopt 7 = 1.0 
for the largest halos (10" - IO^^Mq) and 7 = 1.3 for ha- 
los with masses in the range of 10^^ — IO^^Mq (as favored 
by hydrodynamical simulations which account for the ef- 
fects of baryonic contraction [53] and as observed among 
early- type galaxies |25|). To determine the concentration 
of a halo of a given mass and at a given redshift, c(z, M), 
we use the model of Bullock et al. [55] (alternatively, us- 
ing the mass-concentration relationship of Ref. [37] modi- 
fies our results only slightly) . We also account for halo-to- 
halo variations in concentration, which are taken to fol- 
low a log- normal distribution with A(logj^Q c) = 0.2 [28) . 

Within a given dark matter halo, there are a large 
number of smaller subhalos, each of which contribute to 
the total annihilation rate. In our calculations, we adopt 
two treatments to describe the effect of substructure on 
the annihilation rate. First we use a generalization of the 
analytic model of Kamionkowski et al. [53] , which is cal- 
ibrated to the results of the Via Lactea simulation [30] . 
In particular, we adopt an annihilation rate boost factor 
which follows the radial profile given by their Fig. 4, and 
provides an overall integrated boost of 17 for a Milky 
Way-sized halo [29]. For larger and smaller halos, we 
scale the overall boost factor by M^^^^ [3TJ and adopt a 
radial dependence which follows Ref. [29j, scaled to the 
virial radius of the halo. As a second case, we adopt 
a similar substructure model, but normalized to the re- 
sults of the Aquarius simulation j31, 32 . This substruc- 
ture model yields an overall boost factor that is approx- 
imately 4.5 larger than that of the Via Lactea model. 
Note that these substructure models rely on extrapola- 
tions over many orders of magnitude in halo mass (from 
the current resolution of simulations to the smallest halo 
masses, taken to be 10~^ Mq), and are thus quite uncer- 
tain. On the other hand, because most of the annihila- 
tion boost from substructure occurs in the outer regions 
of halos where the magnetic fields are fairly weak, the 
flux of synchrotron emission is only modestly impacted 
by substructure, in particular when compared to the im- 
pact of substructure on gamma ray emission (neglecting 
substructure entirely reduces the radio flux only by a fac- 
tor of a few, relative to the Via Lactea case) . 

The magnetic fields found in galaxies and and galaxy 
clusters are produced by complicated and not easily mod- 



elled baryonic physics and are difficult to reliably esti- 
mate. They are also rather poorly measured and are 
thus only modestly constrained by observations. In our 
calculations, we adopt a simple and reasonable estimate 
for the distributions of magnetic fields, while acknowledg- 
ing the considerable uncertainties being introduced. We 
will return to this issue in Sec. |IV| For the largest halos 
(lO^'^ — IO^^Mq), corresponding to the scales of galaxy 
clusters, we adopt the following parameterization |33j : 
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where 3/3 (which we take to be 3 x 0.6) is the slope of the 
cluster gas profile at r ^ re and Vc is the core radius of 
the gas distribution (taken to be 5% of the halo's virial 
radius, r^ = 0.05i?vir)- We adopt a central field strength 
of Bq ~ 10 /iG which falls off according to 77 = 0.5. 

For magnetic fields in somewhat smaller, galaxy-scale 
halos {10^'^ — 10^^ Mq) , we start with a double exponential 
model: 



B{r, z) = Bq exp(-i?/i?o) cxp(-z/zo). 



(4) 



with Bq = 7.6^G, Rq = 0.17i?vir, and zq = 0.01i?vir- In 
light of observations which indicate the presence of rather 
strong magnetic fields in the inner region of the Milky 
Way [31], we replace the above parameters within the 
inner several kiloparsecs of galaxy-scale halos by Bq = 35 
^G, Rq = 0.017 i?vir, and zq = 0.006 i?vir (at a given 
location, we adopt the parameter set which yields the 
largest value of B, allowing for a continuous transition 
between the two sets of parameters). 

Although we neglect any isolated (non-substructure) 
halos smaller than 10^'^Mq in our calculations, we note 
that if we had included halos in the mass range of 
101°Mq > M > IO^Mq with 1 iiG fields, they would 
contribute only on the order of 1% of the total radio 
emission. 

In calculating the contribution of dark matter annihi- 
lation products to the isotropic radio background, one 
might worry that magnetic fields in dark matter halos 
may evolve with redshift, potentially contributing less 
from more distant, and thus older, halos. As it turns out, 
however, relatively nearby halos overwhelmingly domi- 
nate this signal (halos at z < 1 provide approximately 
98% of the total contribution to the radio background 
in most of the models considered here). The primary 
reason for this can be seen in the left frames of Fig. [2] 
where we plot the fraction of energy in electron/positron 
dark matter annihilation products which is emitted as 
synchrotron (as opposed to inverse Compton scattering 
with the cosmic microwave background), as a function of 
redshift. Due to the rapid evolution of the CMB energy 
density with redshift, pcMB oc (1-l-z)'*, synchrotron losses 
play only a minor role for energetic electrons in cosmolog- 
ically distant halos. Furthermore, evidence also suggests 
that magnetic fields in high redshift galaxies and galaxy 
clusters are not significantly lower than those found in 
local halos |35j. 
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FIG. 2: The fraction of energy in electron/positron dark matter annihilation products which is lost to synchrotron emission 
(left) and the power-averaged magnetic field strength which contributes to that synchrotron emission (right). In the upper and 
lower frames, we show result for our models for the largest dark matter halos {M > lO^'^i\f0) and somewhat smaller halos 
(lO^'' — 10^'^ Af©), respectively. In each frame, the solid lines denote the results for our default magnetic field and dark matter 
distribution models. The dashed lines represent the range of values found using the variety of models described in Sec. |IV| In 
each case shown, we use the Via Lactea substructure model. 



The diffuse and isotropic spectrum of synchrotron 
emission from dark matter annihilations taking place 
throughout the universe is described by: 
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terested in, energy losses are dominated by inverse Comp- 
ton and synchrotron processes: 
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where pdm is the average cosmological density of dark 
matter, h{z) — ^/Ha + ^m{^ + z)^, and A'^{z) is the 
average squared ovcrdcnsity of dark matter, vifhich is 
determined from a combination of the halo mass func- 
tion, halo profiles, and contributions from substructure. 
dNgyn/dEsyn is the spectrum of synchrotron emission per 
dark matter annihilation, which depends on the spectrum 
of electrons/positrons injected, and on the magnetic field 
strength and other characteristics of the environment in 
which the annihilations take place. The steady-state dis- 
tribution of cosmic ray electrons and positrons can be 
approximately determined by solving the diffusion-loss 
equation: 

= V-[K{x,E,)V^ix,E,)] 

+ JL[t{x,E,)^^{x,E,)]+Q{x,E,), (6) 

where K(x,E^) is the diffusion constant and b{x,Eg^) is 
the energy loss rate. In the energy range we are most in- 
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where p,ad is the energy density of radiation and B is 
the magnetic field strength. Note that in magnetic fields 
stronger than 3.25 /xG, energy emitted as synchrotron will 
exceed inverse Compton losses to the cosmic microwave 
background (at z = 0). 

In large halos, such as those of galaxy clusters, diffu- 
sion can safely be neglected (variations in the magnetic 
field strength and dark matter density occur on length 
scales much larger than electrons and positrons travel 
before losing their energy). In this case, the steady-state 
equilibrium number density of electrons/positrons at a 
given location is given by: 
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where b{Ee) again is the energy loss rate and dNejnj/dE 



is the spectrum of electrons/positrons injected per dark 
matter annihilation. In smaller halos, electrons and 
positrons can potentially diffuse into regions with non- 
negiigibly differing magnetic fields before losing most of 
their energy, which can alter the resulting synchrotron 
spectrum, and the fractions of energy lost to synchrotron 
and Inverse Compton. To estimate the scale of this ef- 
fect, we calculate the typical distance that an electron 
diffuses before losing the majority of its energy to syn- 
chrotron. In a magnetic field of strength B, an electron 
of energy Eg^ will lose half of its energy in a time equal to 
ii„,, = Ee/b{Ee) w 4 X IQi^s X (10GeV/^e)(10AiG/B)2. 
During that time, a typical electron will diffuse a dis- 
tance of Loif ~ \/ K {Ee)tioss{Ee) , which for a diffu- 
sion coefficient of K{Eg) sa 5 x 10^* cm^/s yields Lj^\i w 
1.45kpcx(10GeV/Se)°-^(10MG/B). As such a distance 
scale is not entirely irrelevant in Galaxy-scale halos (but 
is irrelevant on cluster-scales), we will attempt to more 
quantitatively assess the impact of diffusion on the over- 
all synchrotron spectrum and total power that results 
from dark matter annihilations in Galaxy-scale halos. To 
accomplish this, we employed a Monte Carlo technique 
to propagate electrons and positrons generated in dark 
matter annihilation over the region through which they 
diffuse in one energy-loss time (calculated using the com- 
bined magnetic field and CMB energy densities at the 
location at which the annihilation takes place), and then 
obtain the average ratio of the magnetic field to GMB en- 
ergy densities in this region. Employing a diffusion con- 
stant of 5 X 10^* cm-^/s, using both the magnetic field and 
dark matter density models described above, we find that 
10 GeV electrons generated in dark matter annihilations 
deposit approximately 35% less power to synchrotron ra- 
diation than is found when diffusion is not considered. 
For our default parameters, we find that approximately 
56% (86%) of the overall radio signal comes from halos 
larger than 10^^ Mq when using the Via Lactea (Aquar- 
ius) substructure model. For such large halos, the effects 
of diffusion are negligible. We thus estimate the overall 
suppression of the dark matter contribution to the ra- 
dio background due to diffusion to be on the order of 
~5-15%, which is a small effect when compared to the 
significant uncertainties in the dark matter distribution 
and magnetic field models being used. Diffusion will also 
lead to a slight softening of the resulting synchrotron 
spectrum, although this effect is small when compared 
to the uncertainties in the magnetic fields. 

From the steady-state spectrum of electrons and 
positrons, dNf./dEe{Ef,), we can calculate the spectrum 
of synchrotron emission per annihilation: 
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Note that the quantity dN^/dEe is proportional to the 
annihilation cross section and the square of the number 
density of dark matter particles, which cancel with those 
quantities in the first portion of this expression. The 



spectral shape of synchrotron emission is described by 
the function F{x) = a;/^ dzK^/^{z), where i^s/a is the 
modified Bessel function of the second kind. At high and 
low energies/frequencies, the limiting behavior of F[x) is 
given by: 
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We know of no simple closed form to describe the be- 
havior of F{x) at intermediate energies. The critical 
frequency, Vc (and thus the critical energy, E^ = hi'c)^ 
of synchrotron emission is determined by the electron 
energy and by the strength of the magnetic field, v^ — 
BEl X (3ce/47rm3c6). 

With this formalism, we turn in the following section 
to the resulting spectrum and flux of the isotropic ra- 
dio background generated from annihilating dark matter 
particles. 



III. RESULTS 

In Fig. [3J we plot the spectral contribution of dark 
matter to the isotropic radio background for a number of 
dominant annihilation channels and dark matter masses. 
In each case, we have chosen the annihilation cross sec- 
tion to fit the overall normalization of the observed excess 
isotropic background, and have adopted the Via Lactea- 
motivated substructure model (see Sec.pl]). In the upper 
left frame, we see that dark matter annihilations to bh 
predict a spectrum of synchrotron emission that is too 
soft to account for the observed radio background (we 
find fits of x^ = 44.2 and 34.7 for todm =10 and 100 
GcV, respectively, over 12 degrees of freedom). We find 
similarly poor fits for annihilations to other hadronic final 
states, or to gauge or higgs bosons. Although stronger 
magnetic fields could potentially harden this spectrum 
and improve this fit to some extent, no plausible mag- 
netic field models can bring this case into concordance 
with the observed radio spectrum. 

We are therefore forced to conclude that in order for 
dark matter to account for the observed radio spectrum, 
the annihilations must proceed largely to leptons, as 
shown in the other three frames of Fig. |3] In these cases, 
it is much easier to accommodate the observed spectral 
shape. In particular, 10 GeV dark matter particles anni- 
hilating to electrons and 100 GeV particles annihilating 
to taus, muons, or electrons, each provide good fits to the 
observed spectrum {x^ — 12.3, 14.4, 12.1, and 13. 1, re- 
spectively). Several of these cases, however, also predict 
gamma ray fluxes which exceed the isotropic gamma ray 
background, as reported by the Fermi Collaboration |36| 
(see also, Ref. [37] ). In Fig. |4J we plot the spectrum of 
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FIG. 3: The contribution to the isotropic radio background from annihilating dark matter particles with masses of 10 or 100 
GeV, and annihilating to bb, t^t~ , n'^fi~, or e'^e~ . In each case, we have used our default magnetic field models and dark 
matter distributions, and substructure as favored by the Via Lactea simulation (see Sec. |ll|. In each case, the annihilation 
cross section was selected to best fit the observed radio excess. In the case of annihilations to bb, we find a very poor fit to the 
observed spectrum, corresponding to x^ = 44.2 and 34.7 for 10 and 100 GeV masses, respectively, over 12 degrees of freedom. 
We find similarly poor fits for annihilations to other hadronic final states, and for annihilations to gauge or higgs bosons. Other 
choices of the magnetic field and substructure models do not significantly alter this conclusion. In contrast, annihilations to 
leptonic final states can provide a good fit to the spectrum of the observed excess. Annihilations of a 10 GeV dark matter 
particle to taus, muons, or electrons yield a x^ of 36.0, 24.8, and 12.3, respectively. For a 100 GeV mass, these channels yield a 
X^ of 14.4, 12.1, and 13.1. The quality of these fits can change somewhat with variations of the magnetic field model, although 
realistic models universally favor leptonic annihilation channels. 



isotropic gamma ray emission predicted for the same dark 
matter models shown in Fig.lSJ In each case, two lines are 
shown: upper curves corresponds to Bullock-like concen- 
trations, extrapolated down to 10~^ Mq halos, while the 
lower curves use a flat concentration-mass relationship 
below 10^ Mq. Dark matter particles with ~100 GeV 
masses consistently exceed the observed gamma ray spec- 
trum when the annihilation cross section is normalized to 
produce the observed radio background, regardless of the 
employed annihilation channel. Only lighter dark matter 
particles can match the intensity of the observed radio 
flux without violating these constraints.^ 



We note that under the same assumptions pertaining to magnetic 



In Fig. [5j we show the range of dark matter masses and 
annihilation cross sections that can produce the observed 
isotropic radio emission for several annihilation chan- 
nels, and compare this to the relevant gamma ray con- 
straints. In particular, we show constraints derived from 
the isotropic gamma ray background (dashed lines) |36j 
and from dwarf spheroidal galaxies (dotted lines) [38] . 
For the isotropic gamma ray background constraint, we 
require only that the dark matter contribution does not 



fields and dark matter distributions, our code produces results 
that are in good agreement with those of Refs. |13II14| . In partic- 
ular, for a uniform 10 /iG field and NFW profiles with spatially 
uniform boost factors, we calculate radio spectra (and quality of 
fits) which are very similar to those presented in Ref. [l3j. 
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FIG. 4: The contribution to the isotropic gamma-ray background from annihilating dark matter, in the same models as shown 
in Fig. ^ In each case, the upper curve use an extrapolatation to determine the concentration of halos lighter than 10^ Mq, 
whereas the lower curve fixes the concentration below this mass. When these contributions are compared to the isotropic 
gamma-ray background observed by Fermi, we see that annihilations of heavy (ttidm ^ 10 GeV) dark matter particles, or dark 
matter particles which annihilate primarily to hadronic or t^t~ channels, are unable to produce the observed radio background 
without also exceeding the observed gamma ray spectrum. 



exceed the reported error bar in any given energy bin, 
and adopt a flat concentration-mass relationship below 
10^ Mq. This leads to a result that is very similar to 
the "conservative" constraints quoted by the Fermi Col- 
laboration [39]. For the constraints from dwarf galaxies, 
we have estimated the electron channel constraints based 
on those from the muon channel, which have a similar 
spectral shape (the Fermi collaboration does not provide 
results for the e+e~ case). In the case of dark mat- 
ter which annihilates equally to e'^e", /i '*'/.*" and t'^t^ 
(democratic leptons), we estimate the dwarf constraint 
simply by scaling the tau channel constraint by a factor 
of 1/3. In each case considered, the constraint from dwarf 
galaxies is more restrictive than that from the isotropic 
gamma ray background. 

From the left frames of Fig. [5j we see that if the Via 
Lactea substructure model is adopted, gamma ray con- 
straints from dwarf galaxies exclude the vast majority of 
dark matter parameter space which could potentially ex- 
plain the observed isotropic radio excess (only a small re- 
gion near ttidm ^ 5 GeV annihilating to e^e~ survives). 



If we instead consider the Aquarius-motivated substruc- 
ture model, more possibilities appear. In particular, an- 
nihilations to electrons, muons, or to democratic leptons 
can each potentially produce the observed radio back- 
ground without violating any gamma ray constraints. It 
is interesting to note that the range of annihilation cross 
sections required to generate the radio background in the 
e~'"e~ and democratic leptons cases is very similar to that 
predicted for a similar thermal relic {av ~ 3 x 10^^^ 
cm'^/s). 



We emphasize that uncertainties in the magnetic fields 
and dark matter distributions can shift the annihilation 
cross sections required to generate the observed radio ex- 
cess by a factor of a few in either direction. In the next 
section, we will consider the impact of these uncertainties 
more closely. 



23 




Via LacLea Sub. 


^ 


24 




--><*- 


/"^x^ 


25 






■" 


26 






- 



10 ao 

mc« (GeV) 



50 100 



2n 


' 


' 


' 


' 




T T 


Via La 


tea Sub. 


XX 


S4 








/X 


25 










2b 








- 




-33 


Der 


no 


=rati 


-: Leptons 






/" 




Via 


Le 


ctea 


Sub. 




x^ 




-24 








^' 


;f? 


^^/^ 




-25 






^ 


-26 














- 




-23 


T^T 


Aqua, 




- 


-24 








//- 


-25 






.--- " 


'-''.' 








-26 







-23 


M 


'^M 


Aquarius Sub, 


- 


-24 






^^_^ 


'y'^^ 


-25 






-■--■-- X>^ 


- 


-26 











10 20 50 100 

m™ (GeV) 



Democratic I.eptc 
Aquarius Sub. 




10 ao 

m„„ (GeV) 



50 100 



10 20 

m„ (GeV) 



50 100 



10 20 



FIG. 5: The range of dark matter masses and annihilation cross sections which can provide a good fit to the spectrum of the 
isotropic radio excess, using our default magnetic field and dark matter distribution models (solid line, 95% confidence regions). 
No models annihilating to quarks or gauge bosons were found to provide a good fit to the observed spectrum. Also shown are 
gamma ray constraints derived from the isotropic gamma ray background (dashed lines) and dwarf spheroidal galaxies (dotted 
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FIG. 6: As in Fig. |5] but over a range of magnetic field and dark matter distribution models. See Sec.|IV|for more details. 



IV. UNCERTAINTIES 



In the calculations described in the previous sections 
of this paper, we adopted specific but reasonable mod- 
els to describe the distributions of magnetic fields and 
dark matter. Under those assumptions, we showed that 
annihilating dark matter particles could potentially ac- 
count for the observed excess of isotropic radio emission, 
although only for dark matter models within a relatively 
narrow range of parameter space. In particular, under 
those assumptions, the radio background could originate 
from dark matter annihilations only if the dark matter 



annihilates mostly to electrons or muons, and with a mass 
in the range of approximately 4-20 GeV or 14-40 GeV, re- 
spectively. In this section, we consider the uncertainties 
involved in these assumptions and assess whether other 
reasonable models for the magnetic fields and dark mat- 
ter distributions might open any other significant regions 
of dark matter parameter space capable of generating the 
observed radio background. 

In Sec. [n] we described our default magnetic field 
model for cluster-scale and galaxy-scale halos. In partic- 
ular, for very massive halos (M > 10^^ -^0) we adopted 
a magnetic field model described by the parameters 
Bq = 10 /.tG and rj = 0.5 (see Eq. Isl). In this section. 
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we will instead consider the following range for these pa- 
rameters: i?o = 5 — 20 /iG and ry = 0.5 — 0.8. For galaxy- 
scale halos (10^° — 10^'^ M0), we consider the range of 
Bf) = 5.9 — 9.5 /iG, and include models both with and 
without a strongly increasing field in the inner region of 
such halos (see Eq. H and text that follows). We also 
consider halo profiles with inner slopes of 1.0 — 1.2 and 
1.0 — 1.4 for cluster-scale and galaxy-scale halos, respec- 
tively. 

In Fig. [6J we show how relaxing these assumptions im- 
pacts our results. This figure is similar to Fig. [5] but 
here we have marginalized over the magnetic field mod- 
els and dark matter distributions described in the pre- 
vious paragraph. As expected, this leads to broader re- 
gions of dark matter parameter space which could po- 
tentially account for the observed radio excess. These 
variations in our underlying assumptions, however, do 
not qualitatively change our conclusions. Annihilations 
to T~^T~, for example, still cannot provide a good fit to 
the observed excess without also exceeding gamma-ray 
constraints (annihilations to quarks or gauge bosons also 
continue to fail to provide a good fit to the observed radio 
spectrum). Instead, we find that only dark matter par- 
ticles which annihilate significantly to electrons and/or 
muons can potentially account for this signal, although 
with a slightly larger range of masses than were found 
in the previous section. The range of annihilation cross 
sections needed to normalize the synchrotron signal to ac- 
count for the observed radio excess is also now expanded, 
covering av w (0.4 — 20) x 10~^^ cm'^/s for annihilations 
to e+e^, CTw w (3 — 30) x 10^^^ cm'^/s for annihilations 
to fJ.'^n'^ , and CTU ~ (1 — 4.5) x 10"^^ cm'^/s for annihila- 
tions to democratic lepton final states (after accounting 
for gamma ray constraints). Note that in Fig. [6J we do 
not show the constraints from the isotropic gamma ray 
background, as these depend on the halo profile shapes 
which we are now marginalizing over. 



V. DISCUSSION AND CONCLUSIONS 

A number of observations at frequencies between 20 
MHz to 10 GHz T, T'^ have revealed the presence of a 
significant isotropic radio background I2j. While the ori- 
gin of this emission is currently unknown, it lies in signif- 
icant excess of expectations from both Galactic sources 
as well as unresolved extragalactic radio source popula- 
tions. While a variety of astrophysical sources have been 
considered (quasars, supernovae, etc.), none appear to be 
able to account for the observed emission. As previously 
pointed out by Fornengo et al. [131 E] , halos of annihi- 
lating dark matter provide a collection of very faint and 
numerous radio sources, and thus represent an interest- 
ing possibility for the origin of the excess isotropic radio 
emission. 

The spectral shape of the observed radio background 



is quite hard, which limits the range of dark matter mod- 
els which could potentially account for this emission. In 
particular, we find that only annihilations to leptonic 
final states can produce the observed spectrum. Fur- 
thermore, dark matter particles which are heavier than 
^50 GeV consistently exceed gamma ray constraints if 
their annihilation cross section is normalized to generate 
the radio background. After such constraints are taken 
into account, we find that dark matter particles with a 
mass of approximately 4-40 GeV annihilating to e+e", 
12-50 GeV annihilating to n'^fJ.^, or 8-18 GeV annihi- 
lating equally to e~^e~, fJ.'^fJ-~ and t^t~ can provide a 
good fit to the observed isotropic radio excess, depending 
on the magnetic fields and dark matter distributions that 
are assumed. For such models, the annihilation cross sec- 
tion required to normalize the synchrotron signal to the 
observed excess is av w (0.4 — 30) x lO^'^® cm'^/s, which 
is similar to the value predicted for a simple thermal relic 
(crw « 3 X 10^26 cm^/s). 

We point out that if dark matter annihilations are re- 
sponsible for the observed excess radio emission, then 
a significant fraction of the isotropic gamma ray back- 
ground observed by Fermi must also result from dark 
matter. This result is presently consistent with Fermi 
observations of the anisotropy power spectrum for extra- 
galactic dark matter annihilation, which allow up to 83% 
of the extragalactic isotropic background to result from 
dark matter annihilations (in contrast, the contribution 
from blazars is constrained to make up no more than 19% 
of the diffuse gamma ray emission) [30] ■ 

The characteristics of a dark matter particle required 
to generate the observed isotropic radio excess are very 
similar to those required to account for the gamma ray 
emission observed from the Galactic Center [JT], and 
for the synchrotron emission from the Milky Way's ra- 
dio filaments [TH] and the Inner Galaxy (the "WMAP 
Haze") [33j. In particular, each of these signals can be 
explained by an approximately 10 GeV dark matter par- 
ticle which annihilates to electrons, muons, and taus with 
a total cross section of av ^ (1 — 2) x 10~^^ cm'^/s (assum- 
ing a local dark matter density of 0.3 GeV/cm ). This 
choice of mass, annihilation channels, and annihilation 
cross are well matched to those required to generate the 
excess isotopic radio emission (see the lower right frame 
of Fig. pi). This range of dark matter masses is also sim- 
ilar to that required to generate the anomalous signals 
reported by the DAMA/LIBRA [43 , CoGeNT [44] and 
CRESST-H iS. collaborations. 
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